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ABSTRACT 

We calculate the properties of gamma-ray burst spectral lines formed by 
resonant cyclotron scattering in a radiation-driven outflow. Most previous 
models of line formation in gamma-ray bursts are appropriate at the polar 
cap of a neutron star located no further than several hundred parsecs away. 
However, the BATSE brightness and sky distributions indicate that, if bursters 
are galactic, they are located in a corona at distances greater than 100 kpc. 
At these distances the burst luminosity exceeds the Eddington luminosity and 
the plasma in a polar line-forming region is ejected along the field lines. The 
variation of the magnetic field strength and plasma velocity with altitude in 
such an outflow would seem to prevent the formation of narrow features like 
those observed by Ginga and other instruments. However, this is not the case 
because the majority of scatters occur close to the stellar surface, at altitudes 
below z ~ r\, where r\ ~ 10 5 cm is the size of the photon source. Consequently, 
the interpretation of the observed features as cyclotron lines does not rule out 
burst sources in a galactic corona. 

The outflow model predicts both absorption-like and single-peaked emission- 
like features. The latter do not occur in models with static line-forming 
regions, but have been observed in gamma-ray bursts by the Konus and BATSE 
instruments. 
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1. Introduction 



The source of gamma-ray bursts is one of the most hotly debated questions in 
astrophysics ([Lamb 1995t Paczynski 1995| ). These short, intense bursts of high energy 



radiation were first observed by the Vela 3 and 4 satellites on July 2, 1967 (|Klebesadel 



Strong, fc Olson 1973| ; [Bonnel! fc Klebesadel 1996]) . Today, the main question continues to 
be the distance to the sources - are they galactic or cosmological? Interest in this mystery 
has intensified in recent months due to the hard gamma-ray repeater candidate of October 
27-29 1996 and the spectacular follow-up observations to gamma-ray burst GB970228. 
The October 1996 observations are a sequence of four events over a two day period with 
positions in the sky that are consistent with a single source (|Mccgan et al. 1996Q . They 
are the best evidence to date that burst sources can repeat, and that therefore bursting 
does not necessarily lead to the destruction of the source. X-ray and optical observations 
of the GB970228 field offer the best candidates to date for burst counterparts in other 
wavelengths. These observations reveal a transient x-ray source flCosta et al. 1997] ), a 
transient optical point source QGroot et al. 1997a]) , and an extended optical source ( Groot 



et al. 1997b ; Metzger et al., 1997a,b; Sahu et al. 1997 ) whose positions in the sky are 
consistent with that of the burst and with each other. 



One of the strongest pieces of evidence in favor of the galactic point of view is the 
observation of statistically significant absorption- and emission-like features in the spectra 
of some bursts, and the interpretation of these features as cyclotron lines (Mazets et al. 
1981, 1982; |Lamb 1982| ; [Bussard fc Lamb 1982j ; |Mazets et al. 1983| ; [Hueter 1984] |Hueter 
1987|; |Murakami et al. 1988j; |Fenimore et al. 1988|; |Harding fc Preece 1989|; [Wang et~aT 



1989| ; [Alexander fc Meszaros 1989| ; [Lamb et al. 1989] ; |Graziani et al. 1992| ; [Yoshida et al. 



19921 ; IMshimura fc Ebisuzaki 19921 [Barat 1993| ; |Mshimura 1994] ; |Mazets et al. 1996| ; |Briggs| 



1996| ; Briggs et al. 1996|) . Cyclotron lines are formed by the resonant scattering of photons 
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by electrons in a magnetic field, B. The motion of free electrons perpendicular to the field 
is quantized into energy levels ([Landau 1930| ; [Ternov, Begrov, fc Zhukovskii 1966[) , 



E e = Vl + 2nb m e c 2 (1) 

where n is an integer > 0, b = B/B c , B c = 4.414 x 10 13 gauss is the critical field where 
the cyclotron energy Eb = b m e c 2 = 11.6 B u is equal to the electron rest energy, and 
B\2 = B/10 12 gauss. When b <C 1, the energy level spacing is approximately harmonic, i.e. 
E e « m e c 2 + nEs- Photons which possess the right energy to excite an electron to a higher 
Landau level can scatter resonantly into or out of the line of sight, producing spectral 
features. 

Although the observation of spectral features in gamma-ray bursts and their 
interpretation continue to be controversial, the existence of Landau levels and the formation 
of cyclotron lines in the spectra of persistent astrophysical sources with strong magnetic 
fields are well established experimentally. The quantization of electron energies in metals 
is responsible for the de Haas-van Alphen effect - the oscillation of magnetization with 
magnetic field strength ( |de Haas fc van Alphen 1930|) . Cyclotron lines appear in the x-ray 
spectra of about a dozen accretion-powered pulsars ( |Mihara 1995| |Makishima fc Mihara 



. These sources include Her X-l ( |Trirmper et al. 1978Q , 4U 0115+63 ( |Wheaton et al 



|T979D , 4U 1538-52 flClark et al.~H?9"0D , and A0535+26 flCrove et al. 1995| ). Visvanathan 



and Wickramasinghe (1979) made the first detection of cyclotron lines from a white dwarf 
in infrared observations of the AM Her type binary VV Pupis. Since then cyclotron lines 
have been detected in about a dozen AM Her binaries, including AM Her itself (for a 
Chanmugam 199"2j ). Comparing these observations with theoretical models 



review see 



of line-formation reveals much about the physical conditions in the line forming region, 
including the magnetic field strength, the plasma density, and the temperature. 
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The lines reported in gamma-ray bursts are at energies ~ 20 — 60 keV, corresponding 
to magnetic fields Byi ^2 — 6, consistent with a neutron star source. The Ginga satellite's 
observation of harmonically spaced features in three bursts strongly supports the cyclotron 
line interpretation. The features persist for as long as 9 seconds QGraziani et al. 1992] ). 



If the cyclotron line interpretation is correct, then the magnetic field must remain steady 
at least as long as the line is observed. The requirement of a steady field rules out some 
cosmological models, such as coalescing neutron stars, which require the catastrophic 
annihilation of the source. 

Establishing that the observed features are indeed cyclotron lines requires comparing 
the observations with theoretical models that are appropriate for the physical conditions 
expected at the source. Most previous calculations of line formation in gamma-ray bursts, 
such as those of Wang et al. (1989), model the line-forming region as a static slab of plasma, 
threaded with a uniform magnetic field oriented perpendicular to the surface of the slab. 
Such a model is appropriate, for example, at the magnetic polar cap of a neutron star with 
a dipole field, located less than several hundred parsecs from the earth. Lamb, Wang, and 
Wasserman (1990) show that for larger distances the static polar cap model is not valid; 
the burst luminosity exceeds the Eddington luminosity and the radiation force creates a 
relativistic plasma outflow along the field lines. 

However, in order to explain the burst brightness and sky distributions observed by 
BATSE, it has been suggested that, if the bursters are galactic, they are located in a 
galactic corona at distances of 100-400 kpc (for a review, see [Lamb 1995| ). Two models 



have been suggested that are appropriate for these distances: an equatorial model and an 
outflow model. 



In the equatorial model, the line-forming region is located at the magnetic equator, 
where field lines parallel to the stellar surface confine the plasma magnetically (see, e.g. 
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Lamb 1982| ; |Katz 1982| ; Zheleznyakov and Serber 1994, 1995). Freeman et al. (1996) fit 



both this model and a static polar cap one to the two observed spectra corresponding to 
the time intervals SI and S2 during which lines are observed in GB870303. A joint fit to 
the two intervals, using models with a common magnetic field and column depth, but not 
a common viewing angle, marginally favors the equatorial model over the static polar cap 
one. 

In the outflow model the radiation force on the plasma at the polar cap ejects the 
plasma relativistically along the field lines. Table 1 compares the existing calculations of 
the properties of cyclotron lines formed in such an outflow. The variation of magnetic 
field strength and plasma velocity with altitude would seem to prevent the formation of 
narrow scattering lines. However, Miller et al. (1991, 1992) show that narrow lines can be 
formed in an outflow at the second and third harmonics, which they approximate as due to 
cyclotron absorption. Chernenko and Mitrofanov (1993) calculate the properties of the first 
harmonic line, also approximating it as due to absorption, and find that the formation of 
a narrow line is possible. However, the absorption approximation is not valid for the first 
harmonic; multiple photon scatters must be taken into account. The first calculation to 
do so was Isenberg, Lamb, and Wang (1996) using a Monte Carlo radiative transfer code. 
This calculation includes multiple scattering at the first three harmonics in a plasma with 
electron temperature ksT = 0.25 Eb- It shows lines with equivalent widths comparable to 
those observed by Ginga. 

The emerging photon spectrum is very sensitive to the velocity profile (3f{z) - the 
bulk flow velocity as a function of the altitude z above the stellar surface. An electron with 
velocity (3 along the field is overtaken by photons with [i > (3, where \x = cos 9, and 9 is 
the angle between the photon's direction of travel and the magnetic field. In the electron's 
rest frame, these photons have orientations // > 0; scattering with them accelerates the 
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electron. Similarly, photons with jj < (5 have rest frame orientations /i r < and decelerate 
the electrons. Consequently, the three previous calculations assume that the flow velocity at 
a given altitude is approximately the mean value of the component of the photon velocity 
along the field; that is 

SQ(ji,z)dfi 

where Q(n, z) is the angular distribution of photons at altitude z. This profile is a good 
approximation to the one Tsygan (1981) and Mitrofanov and Tsygan (1982; hereafter MT) 
derive from the radiation force F ra d(z). Their derivation includes scattering at the first 
harmonic only, with k^T = 0. It is based on the unscattered photon distribution; i.e., the 
calculation of the force at z does not take into account the reduction in photon flux near 
the line center or the angular redistribution of photons due to scattering at lower altitudes. 
This approximation is valid at small optical depths, but not at the depths required to 
reproduce the observed cyclotron lines. Miller et al. (1991, 1992) and Isenberg, Lamb, and 
Wang (1996) also include processes in their radiation transfer calculations, such as higher 
harmonic scattering and finite temperature, which are not included in MT's calculation of 
F ra d. Thus, in the earlier models, the hydrodynamics is not consistent with the radiation 
transfer. 

In the current work, we present the first self- consistent calculations of the properties 
of cyclotron lines formed in an outflow. In addition to the physical processes included in 
Isenberg, Lamb, and Wang (1996), we add continuum scattering, finite natural line width, 
and photon polarization. We describe the geometry and the physics of the line forming 
region in sections |2| and ^| respectively. The electron velocity profiles are discussed in section 
f|. We present the Monte Carlo spectra in section [| and discuss their implications for 
gamma-ray bursts in section |6|. 
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2. Spatial Geometry of the Line-Forming Region 

We adopt a cylindrical line-forming region, as shown in figure [l|. Photons are injected 
at the base of the cylinder, which could correspond to a hot spot with radius located on 
the surface of a neutron star. The angular size of the hot spot is represented by 



p {z) = cos 6 (z) = ^ 1 + ~ 2 y /2 , (3) 

where z is the altitude above the surface and the tilde indicates quantities measured in 
units of Th- 

We assume the star has a dipole magnetic field with strength B , where the subscript 
o indicates quantities measured at the stellar surface. We also assume that the hot spot 
is located on the stellar surface, at the magnetic pole. The field strength decreases with z 
according to 



where is the stellar radius. The field lines that are rooted in the hot spot form a tube 
with a circular cross section with radius r(z). Conservation of magnetic flux requires that 
r(z) increase with altitude according to 



r - w = ,S: = ( 1 + ir*, (5 ) 



, B(z) V R 

i.e. the field lines flare outwards. For simplicity, we assume the field lines are parallel to 
the cylinder axis at all z; we ignore the flaring in our Monte Carlo simulations, except 
when calculating the plasma density. The characteristic distance for changes in r(z) is 
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~ (df/dz)^ 1 ~ i?*. But, as we shall show, few scatters occur above z ~ <C R*. Therefore, 
the effect of this assumption on the emerging photon spectrum is negligible, since most of 
the scattering takes place close to the surface where r(z) is approximately constant. 

The line-forming region contains a plasma in which the number density of scatterers 
(electrons and positrons) is n e (z). The plasma can consist of ionized electron-proton pairs 
that have been swept up from the stellar surface or electron-positron pairs that have been 
created by the processes 7 — > e~e + and 77 — > e~e + . Our Monte Carlo simulations assume 
the number of scatterers is conserved; i.e. the simulations apply to line-forming regions 
where electron-proton and electron-positron plasmas originate at the surface. The effects 
of pair creation and annihilation above the surface is beyond the scope of our current 
calculations (although, see section || for a general discussion; see also |Vitello fc Dermer 



1991| ; |Coppi fc Lamb 1992| ). Because the electron temperature k^T e <C Eg, the electrons 



are restricted to move along the field lines (see section |3.2| ). As the photons travel through 
the plasma they can scatter with the electrons (and positrons), accelerating them. We 
denote the velocity along the field lines of an individual electron (in units of c) by (3 and 
the bulk flow velocity by (3f{z). 

Photons which cross the lateral surface of the cylinder escape the line-forming region 
and reach the observer. In the absence of scattering, a photon injected at the origin escapes 
when fi (z) = fi. The resulting angular distribution of photons is 



Q(jm,z) = Q(n) rj 



H - n {z) 



(6) 



where Q(/i) is the angular distribution at the stellar surface and rj is the unit step function. 
For simplicity, we use semi-isotropic injection (Q(fj) = 1) throughout the present work. 
Note however that some models of spectrum formation in gamma-ray bursts suggest that 
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the photons are beamed along the field lines (see e.g. |Schmidt 1978| ; Baring fc Harding" 
1993| ; [Harding fc Baring 1994| ; [Harding 1994| ; Section ^ below). This possibility should 
be considered in future calculations of line-formation in an outflow. Also note that the 
distribution given by eq. (§) for photons injected at the origin of the hot spot is also the 
angular distribution along the cylinder axis when photons are injected uniformly across the 
disk. Substituting this distribution into eq. (El) yields the approximate velocity profile 



P F {z) 



1 + Ho(z) 



(7) 



The few photons which reach z > 30 are also permitted to escape without further 
scatters; for ~ 0.1.R*, the plasma is extremely tenuous above this altitude and the 
cyclotron energy less than the energy threshold of most gamma-ray burst detectors. We 
assume that photons which scatter downwards and return to the surface thermalize, i.e. 
they are absorbed by nonresonant inverse magnetic bremsstrahlung. 



3. Physics of the Line-Forming Region 

In this section we describe the basic physics that govern the properties of cyclotron 
lines formed in outflows: the characteristic time scales (section |3.1| ), the one dimensional 
electron momenta distribution f(p) (section |3.2|) , the plasma density profile n e (z) (section 
|3~3|) , the scattering kinematics (section |3]4]), the scattering profile <fi(E,fi, z) (section [3.5|) , 
and the gravitational red shift (section |3T6|) . We discuss the electron velocity profile (3p{z) 
separately, in section |J 



3.1. Characteristic time scales 



Five time scales characterize the physics of the line-forming region: the radiative 
lifetime t ra d, the electron-photon relaxation time scale t e7 , the electron-proton and 
proton-proton relaxation time scales t ep and t pp , and the particle escape time t esc . For 
simplicity, we estimate the first four time scales for a static [y F = (1 — /3f,) -1 / 2 = 1] 
line-forming region. This approximation is suitable near the stellar surface, where < 1.2. 
All these time scales are short compared to the burst duration ~ 1-10 s. 

The radiative lifetime of an electron in the nth Landau level is just 



r 

where 



nnl = A = 3 x 10- 16 s n- x B^ , (8) 



T„ = = 2.6 x 10- 3 keV nB\ 2 (9) 

is the natural line width for the nth harmonic and a is the fine structure constant. 



We determine the time scale for the electrons to come to equilibrium with the photons 
by calculating the Fokker-Planck coefficient 



where 



6ptt(fji- ii s ) E 



(11) 
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is the change in momentum during a scattering, the subscript s denotes parameters of the 
scattered photon, and 

n x (E, n, z) dE^- = Q(/i, z) n x (E, z) dE^ (12) 

271 ATT 

is the density of x-ray photons at altitude z with energy between E and E + dE and 
direction between Q and Q + dQ. Assuming first harmonic scattering dominates the 
radiation transfer and taking the limit of zero natural line width, the polarization averaged 
scattering cross section is, to lowest order in (3 and b, approximately 

W. w 32 aT — 6{E ~ Eb) — — (13) 

(see Section [3.5|) , where ot = (87r/3)[e 2 /(m e c 2 )] 2 = 6.65 x 10~ 25 cm 2 is the Thomson cross 



section. Substituting eqs. ([11]) - fll3|) into eq. ( |10| ) yields 



5p 2 \ 137r m e c 3 , _ . „. 

tTT-^-n^C^fl) . (14) 



5t / 160 a 



The equilibration time scale is 



If the electron column depth is sufficiently low so that photons scatter at most one time 
before escaping the line-forming region (the single scattering approximation), then at all z 
the photon spectrum is approximately the same as the spectrum initially injected at the 
hot spot. For a power law spectrum with spectral index s, 



13 



n x (E,z) =n x {E) = x E~ s 

nric f (s) 

= 6.6 x 10 27 keV cm" 3 J*g- ( , (16) 

rl, 5 £00 VI keW 

where Th,5 = f/i/10 5 cm, L x ^q = L x /10 erg s _1 is the luminosity of x-ray photons with 
E x < E < E 2 , and 



£(s) = E l ~ s dE . (17) 



Throughout this paper we use Ei=l keV and £ , 2 =1000 keV. Substituting eqs. (^) and 
(0) into eq. (|15|) and choosing s = 1 [£(s) » 1000 keV], 

t e7 ~ 5 x 10- 15 s L-% < 5 B 12 . (18) 
By similar arguments, the relaxation time for electron-proton scattering is 

t ep ~ 3 x ID" 6 s n- p \ 7 Tf 2 1 , (19) 

where = fceTg/l keV, the Coulomb logarithm is lnA ep « 16.2 + 2 In (T^-B^) ( |L anger 



1981 ), and np,i7 = n p /10 17 cm -3 is the proton number density (= n ei n = n e /10 17 cm -3 
in a pure hydrogen plasma). We choose a characteristic density, 10 17 cm~ 3 , where the 
line-forming region is optically thick in the line core but optically thin in the wings (see 
section |37| below, esp. eqs. [41] and [47] and associated discussion). The relaxation time 



for proton-proton scattering is 



i PP ~2x 10- 7 s n-\ 7 Tl /2 P^j*) , (20) 
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where lnA pp « 15.9 + 0.5 ln(T fe n~\ 7 ) ([Spitzer 1962ft . Since the electrons are confined to 
motion along the field lines, they merely exchange momentum when they scatter with each 
other, with no change to the overall momenta distribution. 

Finally, the time required for an electron to escape the region near the surface, where 
most scatters take place, is 



t esc ~- = 3x 10" 6 s . (21) 

c 

Comparing the time scales reveals a great deal about what physical processes 
are important in the line-forming region. Since t ei <C t ep . m radiation processes 
dominate collisional processes. The reason for this is clear: the density of electrons 
is small compared to the density of resonant photons. While the electron number 
density in the outflow model is 10 17 — 10 19 cm -3 , the density of resonant photons is 
« n x (E B ) £]= 4 x 10 23 cm" 3 L xA0 T fc 1/2 r^j, where 



E n d = nE B J = 0.73 n£ 12 T fc 1/2 (22) 
y m e c z 

is the Doppler width associated with the n th harmonic. Since t ePjPP may be ^ t esc , the 
particles may not spend enough time in the line-forming region to reach equilibrium by 
particle-particle collisions. However, t ej <C t esc , so there is ample time for the electron 
momenta distribution to reach equilibrium by electron-photon scatterings, as we discuss in 



more detail in Section 3.2. 



Since particle-particle collisions are rare, radiative processes determine the population 
of excited Landau states. As eqs. (^j) and flTS|) show, for most situations of interest for 
line formation in gamma-ray bursts the radiative decay time of the excited states is shorter 
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than the mean time between electron-photon scatters. Consequently the excited states are 
not radiatively populated. Thus, we assume throughout this work that an electron is in 
the Landau ground state (n = 0) before and after a scattering. Note, however, that this 
assumption breaks down for low field strengths and high luminosities where Bf 2 L~^ ^0.1. 

The low rate of particle-particle collisions suggests that pure absorption (i.e. 
nonresonant inverse magnetic bremsstrahlung) is rare. A more careful analysis confirms 
this. For Eb 3> ksT e , the absorption probability per scattering is the ratio of absorption to 
scattering cross sections 



Pa = 


~cr a ~ 




~cr a ~ 











2ira 



lm P c 



2 h 3 c 3 n T 



B=0 



E 



B 



E% 



1.5 x 1(T 10 n pA7 B\ 7/2 



'12 



(23) 



( [Nelson, Salpeter, fc Wasserman 1993] ) . If the line wings are optically thin (see Section |3.5|) , 
a photon in the line core can escape the line-forming region by scattering to the wings. A 
core-wing transition occurs in ~ 1/a scatters, where 



St- = 1.8 x 1(T 3 B 12 T k l/2 



(24) 



is the dimensionless natural line width ( Wasserman Salpeter 198U ). The total absorption 
probability for a resonant photon is, therefore, a~ 1 P a ~ 10~ 7 n p n B^^ 2 T^ 2 . Since this 
probability is 1 throughout the current work, we do not include pure absorption in our 
calculations. 



-16- 



3.2. Electron Momenta Distribution 

The time scales that govern the relaxation of electrons in a dynamic, galactic corona 
model are dramatically different from those of a static, galactic disk model where the 
photon densities are much lower and the electrons do not escape the line-forming region. 
Consequently, determining the distribution of electron momenta along the field lines in 
an outflow requires some care. Two questions must be answered: what is the equilibrium 
distribution and are the electrons in equilibrium — i.e. do they possess the equilibrium 
distribution? 

In the frame of reference moving with the flow, the equilibrium distribution is 
Maxwellian. We demonstrate this by showing that the electrons and photons constitute 
a canonical ensemble in equilibrium. In this system the electrons compose a subsystem 
in thermal contact with a heat bath made up of the photons. To constitute a canonical 
ensemble, the particles must satisfy two requirements. First, the number of particles in the 
subsystem must be much smaller than the number of particles in the heat bath. Second, 
the sum of the energies of the system and the heat bath must be constant (see, e.g., [Landau 



k, Lifshitz 1980Q . Since n e <C n x E\ and energy is conserved in electron-photon scattering, 



both of these conditions are satisfied. Thus, the probability that an electron has an energy 
(p 2 + m^c 4 ) 1 / 2 is a Maxwell-Boltzmann distribution. 

The electrons relax to equilibrium by electron-photon scattering in a time ~ t e7 . 
Note that there is no requirement that the photon distribution be thermal or that there 
be significant electron-proton scattering for the electrons to reach equilibrium. Since 
t ei <C t esc , the electrons in the outflow have an equilibrium distribution, except in the tail 
of the Maxwellian; the electrons do not remain in the line-forming region long enough to 
scatter a sufficient number of times to populate the tail. Since the number of electrons in 
the tail is small, however, they do not have a significant impact on the properties of the 
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emerging radiation. Therefore, in our calculations we assume an equilibrium distribution at 
all electron momenta. 

The electron temperature at equilibrium (i.e. the Compton temperature) Tc is equal to 
the temperature of the photons T 7 = (dS '/ dU)' 1 , where S and U are the entropy and energy 
of the photons, respectively. The Compton temperature can be determined microphysically 
by summing the energy transferred from the radiation to the plasma in each scatter; at the 
Compton temperature this sum is zero. Using the single scattering approximation, Lamb, 
Wang, and Wasserman (1990) calculate Tc in the limits 6,/3 < 1 for a static plasma and 
first harmonic scattering only. 

We adapt their model to a dynamic plasma by finding the Compton temperature T c , 
where the superscript f indicates quantities measured in the frame of reference moving with 
the flow. We then boost the momenta distribution back to the lab frame. From eq. (^), eq. 
(|I6|), and the Lorentz invariance of E~ 2 n x (see, e.g. [Landau fc Lifshitz 1975| ), 



<(£V) = 7 ^ , (25) 



7f(1 + (3 F ^) 



Substituting ( plf) into the formula of Lamb, Wang, and Wasserman (1990), 



. T f = Eb Jg (2 + 7// 2 + 5// 4 ) ni(E B , //) d^ 
B c 10 J+l [1 + (2 + s)^ + (s- 3)/i f4 ] nl(E B , /i f ) d^ ' 1 ' 

and using eq. ([/]) for the velocity profile gives ksT c w 0.25Eb for all z. We use this value 
of the Compton temperature throughout the present work. For a Maxwellian distribution 
with temperature /cgT f •jp 2 m e c 2 in the frame moving with the flow, the lab frame 
distribution is approximately Maxwellian, with the peak aX p = pf = If Pf and 
temperature T e = 7^T e f . The approximation is a good one near the peak, but breaks down 
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in the tail. However, as we point out above, the shape of the distribution in the tail does 
not have a significant impact on the properties of the emerging radiation. 

Monte Carlo calculations by Lamb, Wang and Wasserman (1990) and Isenberg, 
Lamb, and Wang (1997) show that when higher harmonics, multiple scattering, and the 
geometry of the line-forming region are included, Tq varies from the analytic value (|26|) 
but is always between 0.2Eg and 0.7Eb- For the line energies observed by Ginga, the 
full-width half-maximum of a line formed in an outflow is relatively insensitive to the 
temperature since the contribution of thermal broadening to the the line width is small. 
This contribution is just the Doppler width E\ = 3.2 keV for Bu = 1.7, T| = 5.0, and 
(3f = 0.5. By comparison, the broadening of the line due to the variation in the field with 
altitude is AE B = \dE B /dz\ = 6.0 keV for the same field strength. Therefore the total 
line width is « [(AE B ) 2 + (E}) 2 ] 1/2 = 6.8 keV w AE B . We show in Section £5| that the 
optical depth, and therefore the equivalent width, is also insensitive to the temperature. 
Therefore, eq. (p6|) is sufficiently accurate for determining the properties of the emerging 
spectrum. 



3.3. Plasma Density Profile 

Conservation of matter in an outflow requires that in the absence of plasma sources or 
sinks, n e (z) f3p{z) r 2 (z) is a constant. The conservation of magnetic flux determines the 
radius of the outflow r(z), according to eq. @. Therefore, if the plasma at the surface has 
density n eo , the density at altitude z is ( Miller et al. 1991] ) 



n(z)=n lE^ElA (27) 
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3.4. Scattering Kinematics 

A resonant scatter may be thought of as an absorption of a photon by an electron, 
which excites the electron to a higher Landau level, followed by the emission of one or more 
photons. The absorption conserves energy and momentum along the magnetic field; the 
portion of the photon's momentum perpendicular to the field that is not transferred to the 
electron is absorbed by the field. Combining the conservation laws for a transition from 
the to nth Landau state, and ignoring the finite natural line width, yield the resonant 
condition (c = m e = 1 throughout this section) 



1 + E r - ^Jl + 2nb + fi r2 E r2 = , (28) 

where the superscript r indicates quantities measured in the frame of reference where the 
electron is at rest prior to scattering. Solving eq. for the rest frame resonant photon 
energy yields 

= 2 »" . (29 ) 

l + y/l + 2nb (1 - u r2 ) 

Transforming to the lab frame (the frame of reference where the magnetic dipole is at rest), 

En = ( f : , ; ■ (30) 
7(1- Pfj.) 

We plot the first harmonic lab frame resonant energy as a function of electron 
momentum p in figure ^|. As the figure shows, for a given photon energy and < 1 there 
are two roots for the electron momentum; p_ and p + denote the lower and upper roots 
respectively. Physical solutions to the resonant condition do not exist above the cutoff 
energy 
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VI + 2nb-l . . 

E ° = n—n (31) 



( Paugherty fc Ventura 1978 ; [Harding fc Paugherty 1991 ; Wang, Wasscrman, fc Lamb 



|1993 ). To understand the physical origin of the two momentum roots and the cutoff energy, 



consider a frame of reference (indicated with a prime) which is moving with a velocity (3 = p 
with respect to the lab frame. In the primed frame, the photon is moving perpendicular to 
the magnetic field (// = 0), since 

✓ = (32) 
1-/1(3 

In this frame the resonant condition is symmetric with respect to reflection in the plane 
perpendicular to the field. Thus there must be two equal and opposite solutions for p'; 
i.e. the E' n vs. p' curve is symmetric about p' = and p' + = —p'_, as shown in figure 
For non-zero electron momentum the photon appears redshifted in the electron rest frame. 
But at p' = 0, fi r = p! = 0, there is no redshift, and the resonant energy has its maximum 



value E' n = E' c = E r n = a/1 + 2n6 — 1. Boosting p' ± and E' c back to the lab frame gives the 
momentum roots and cutoff energy when p ^ 0. 

Figure |2] illustrates that the scattering kinematics in an outflow are noticeably different 
from those in a static plasma. For a given p, the line center energy E ctr is equal to the 
resonant energy corresponding to (3 = (3p- In a static plasma (/3p = 0), E ctr w for all 
/i. The cutoff energy gives the scattering cross section a strong asymmetry near p = 
( [Lamb et al. 1989|) but plays little role at larger p where the cutoff energy is far above the 



line center. At these angles, the p_ root dominates the scattering. However in an outflow 
with, for example, (3p = 1/2, there is a large variation in the line center energy with photon 
direction, as expected from eq. (|30|). The effects of the cutoff energy are the most dramatic 
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at jj, = (3p, where E ctr = E c . p_ dominates the scattering at larger values of \i while p + is 
dominant at smaller ones. 



3.5. Scattering Profile 

In calculating the scattering profiles we assume that the vacuum contribution to the 
dielectric tensor is much larger than the plasma contribution. This is true provided 



w 



4 x KT 5 n e ,i7 £f 2 4 < 



Eb 
~E 



(33) 



where w = (Twu p /Eb) 2 is the plasma frequency parameter, u p is the plasma frequency, and 



5 is the magnetic vacuum polarization parameter (see |Adler 1971| ; [Wang, Wasserman, fc 



Salpcter 198SQ . For the fields and densities considered in the present work, eq. (|33D breaks 
down above z ~ R* where the field strength becomes small. However, the condition is 
satisfied near the stellar surface where most of the scattering takes place. This condition 
also breaks down very near the line center, where E ~ Eb- However, the effect on the 
emerging spectrum is negligible since the line properties are determined primarily by the 
scattering profile in the line wings (see below). 

The scattering cross sections we use are valid in the limit 



m P c 



(34) 



which holds throughout the present work. In this limit the scattering cross section 
approaches the classical magnetic Compton cross section. When evaluated in the 
pre-scattered electron's rest frame and averaged over the azimuthal angle (appropriate for 
azimuthally symmetric line-forming regions such as those studied here) this cross section is 
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given by flCanuto, Lodenquai, fc Ruderman 1971| ; |Herold 1979| ; |Ventura 1979| ; |Wasserm"an 
Salpeter 1980| ; |Wang, Wasserman, fc Salpeter 1988| ; [Lamb, Wang, fc Wasserman 199Q ; 
Harding fc Daugherty 199l| ; Praziani 1993| ; |Isenberg, Lamb, fc Wang 1997] ): 



da r dal 



+ 



do\ 



dn: da: dn 



The first term on the right is the continuum part of the cross section 
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(36) 



The symbols at the right of eq. (^) indicate the polarization mode, parallel or 
perpendicular, of the initial and scattered photon. The resonant part of the cross section is 



d£Z 



9 aTTn e c 2 
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E B (Er + E B y 



X 
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(- 

1 (- 



(37) 



The corresponding polarization-averaged cross sections are calculated by summing eqs. (^6| 
and (|37|) over the final polarization states and averaging over the initial. 



The contribution of n > 1 harmonic scattering to the scattering profile can be 
calculated approximately by treating higher harmonic scattering as the absorption of a 
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photon followed by the emission of one or more photons. Then the total rest frame cross 
section 



oo 

* r ~E<> (38) 

n=0 



where 



o 



(the sums are over the final polarization states), and 



\n-l 



a 



n>l 



3 W<r r fen _ a (n 2 /2)- x 

4 a ! fn-1)! 



r„/2vr 



{E r _ £r)2 + ( r B/2)i 



r2 



(40) 



is the total cross section for absorption at the nth harmonic (Daugherty and Ventura 1977; 
Fenimore et al. 1988). The approximation ([38|) is valid provided the line forming region is 
not optically thick in the line wings of the first harmonic (see below; also [Wasserman fc 
galpeter 1980); [Lamb et al. 1989|; [Isenberg, Lamb, fc Wang 1997|), i.e. 



7i £ 1/a , (41) 

where T\ is the polarization and frequency-averaged first harmonic optical depth for photons 
moving parallel to the magnetic field (see below), and a is given by eq. ([24]). Eq. ( ]4"T| ) holds 
for all cases we study in the present work. 
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From the Lorentz invariance of the optical depth (see, e.g., Rybicki fc Lightman 1979| ), 



the lab frame cross section a n is related to a r n by 



(i - m < 



(42) 



We average the lab frame cross sections over the electron momenta distribution and divide 
by the Thomson cross section to obtain the scattering profile for the nth harmonic 



The line can be divided into the line core (|x n /u| <C 1) and the line wings (|x n /u| ^> 1), 
where x n = (E — E n )/E% is the dimensionless frequency shift (Wasserman and Salpeter 
1980). In the line core, the electron distribution dominates the profile so that n oc f(x^/\fi\). 
In the wings, the tail of the Lorentzian distribution dominates so that (ft n oc ax~ 2 . We refer 
to the wing at energies below the line center as the red wing and the wing at energies above 
the line center as the blue wing. Wasserman and Salpeter (1980) show that for the first 
harmonic, the core-wing boundary appears at \xi/fi\ « 2.62 — 0.191n(100 a/u) . 

The optical depth at the nth harmonic is equal to 



The first harmonic optical depth along the line of sight with the shortest escape path 
characterizes the radiation transfer. In the outflow model this line of sight is along the field, 
due to the change in the cyclotron energy with altitude: a photon can escape the line core 
by diffusing to an altitude where it is in the wings. If the characteristic distance a photon 
climbs before escaping the line core 




(43) 




(44) 
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Az P 



2=0 



dEs/dz 

then the energy averaged first harmonic optical depth along the field is 



(45) 



(46) 



2=0 



Averaging over the photon polarization and adopting the velocity profile from the model of 
MT and a Gaussian electron momentum distribution yields 



S ,i2(l-4.5r fc /i^) 1 ; 

For rh/R* = 0.1, eq. fl47D gives ri ~ 149 n e0j i 7 77^5 S~J 2) an d the condition eq. fl41"|) 
implies n e0j i 7 ^ 4 r^5 X^ 2 . Comparing eq. (f47|) for rj to the corresponding optical depth 



for a static slab line-forming region is instructive. In a slab with a uniform magnetic field 
parallel to the slab normal, the optical depth along the field is T\ = 1500 N e>2 2 B± 2 T k 
( [Lamb, Wang, fc Wasserman 1990|) , where N e ^2 is the column depth in units of 10 22 cm~ 2 . 
The optical depth of an outflow is about an order of magnitude smaller than that of a slab 
with a comparable column depth, i.e., a slab whose column depth is equal to the radial 
column depth n eo rh of the outflow. The optical depth is smaller for the outflow because the 
distance a photon needs to travel to escape the line core is smaller than r^; it can escape by 
climbing to an altitude where the cyclotron energy is less than the photon energy by more 
than a Doppler width. Note that in the outflow, the optical depth is independent of the 
electron temperature, in contrast to the slab where T\ ~ T fc ~ 1/2 . The reason for this is that 
in the outflow the decrease in dr^jdz with increasing temperature is offset by the greater 
distance a photon needs to climb to escape the line core. 
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3.6. Gravitational Red Shift 

A photon that travels from the stellar surface (z = 0) to an observer at z = oo 
experiences a gravitational red shift 

E(0) - E(oo) _GM* / M* \ ( R* y 1 

E(0) c 2 R* ' \M ch ) llO 6 cm/ ' 1 ' 

where M* is the stellar mass and M c h = 2.8 x 10 33 g is the Chandrasekhar mass. Because the 
red shift changes the photon energy by only « 20%, we do not include it in our calculations. 
The red shift reduces the cyclotron line width as well as the energy of individual photons, 
so the red shift narrows the lines. Since we wish to determine whether narrow lines can be 
formed in an outflow, ignoring the red shift is a conservative assumption. 



4. Electron Velocity Profiles 

The net force on a particle in an outflow is the sum of the radiation and gravitational 
forces, F = F ra d + F s . For the luminosities appropriate for a gamma-ray burst in the 
galactic corona, the radiation force is much larger than the gravitational force, so that 
F f» F ra d. Since the electrons are restricted to motion along the field lines, only the 
component of the radiation force parallel to the magnetic field is of interest. In Section [4. 1| 
we use the approximate, analytic model of MT to calculate the radiation force along the 
field line. We then integrate the equations of motion to determine the velocity profile. In 
Section we use a Monte Carlo radiative transfer code to determine the velocity profile 
self-consistently. 
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4.1. Analytic Model of Mitrofanov and Tsygan 

MT calculate the radiation force in an outflow using a single scattering approximation. 
They show that when first harmonic scattering dominates the radiation transfer, and 
r r — > 0, the component of the radiation force along the field line is 



E 2 r 

■nut T^OT /•.'/,- J (// - 3) 



i-/W V(i-W 



r^W nx [whty M ) dfl ' (49) 

where Pi is the degree of linear polarization, i.e. the fraction of perpendicular mode photons 
minus the fraction of parallel mode photons. Because eq. is linear in the scattering 
cross section, the radiation force for the polarization averaged cross section is the same as 
the radiation force for polarized cross sections with Pi = 0. Substituting the power law 
photon spectrum we use in the present work [eqs. (0) and fll6|)1 into the radiation force 
gives 



Frad ~ — TV -rUW 1 ' (50) 

where 



J fl 



1 + 



l-Pu 



+ Pl 



7 2 (l-/3/i) 5 



(51) 



Eq. (|50D is equivalent to 
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rad 



8 



1-s / tt„\ \"! 



where the upper expressions in curly brackets are for an electron-proton plasma and the 
lower expressions are for an electron-positron plasma. 

By integrating eq. (0), we calculate the velocity as a function of altitude for an 
electron in an electron-proton plasma, injected at the origin with (3 = 0. We use a 1/E 
photon spectrum with Pi = 0. We plot the result in Figure |3|. The curve shows three distinct 
regimes which MT call free acceleration, kinematic restriction, and energetic restriction. 
At any altitude z there exists a velocity f3(F=o) where the radiation force is equal to zero. 
Near the surface, (3(f=o) ~ 0.5, which corresponds to 7(f=o) ~ 1.15. In the free acceleration 
regime [3 < (3(f=o) and the electron accelerates rapidly. 

When (3 reaches (3(f=o) the electron enters the regime of kinematic restriction. The 
distance it travels before reaching this regime is 

ZKR ~ dWrfi ~ ( 10 -« / Ukv} ^ ^ooo kev 2 - J ■ (53) 

In the kinematic restriction (3 ~ (3(f=o)- The electron requires a small radiation force to 
overcome the force of gravity and keep pace with the increase of (3(f=o) with altitude. 
However the difference between the actual velocity and the zero force velocity is ~ one part 
in 10 6 . 



As the electron continues to climb, the radiation becomes more diffuse and eventually 
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carries too little energy to maintain the particle at the zero force velocity; i.e. 



d 7(F=o) > dn _ , 54) 
dz dz 

The particle enters the energetic restriction regime and the Lorentz factor levels off to 
7 = Imax- For z ^> 1 and s — 1, eq. flBTD becomes 



h « J- 

1 2F 



1 3 



< i > (55) 



2P 7 



so 7(f=o) ~ (2/3) 1//2 2;. Substituting 7 ~ z and eq. into eq. Q53D, the energetic 

restriction begins at 

J 100 1 rV3 -1/3 
[ 1000 J 

We caution that the discussion of the energetic restriction is highly speculative. Since the 
cyclotron energy in this regime is below burst detector thresholds, we do not know the 
relevant part of the photon spectrum observationally. 

Figure ||] and eqs. ( |53f ) and (j5f|) clearly show that the electron is in the kinematic 
restriction throughout the region of interest for line formation, except for about a millimeter 
of free acceleration near the surface. Consequently, we assume F ra( i = for the remainder 
of the present work. 



Eq. (|52| ) gives the acceleration of an individual electron. For a group of electrons whose 
velocities are distributed thermally about the flow velocity, the radiation force should be 
averaged over the velocity distribution. Setting the average force to zero and solving for 
the flow velocity yields a velocity almost indistinguishable from the zero force velocity 
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of an individual electron, i.e. (3p ~ (3(f=o)- There are two reasons for this. First, the 
corrections to the velocity due to temperature are ~ kBT e /m e c 2 -C P(f=o)- Second, the 
Maxwell-Boltzmann distribution is symmetric, so the corrections due to electrons with 
velocities on each side of the bulk flow velocity tend to cancel. 

Figure ^ shows the zero force velocity profile given by eq. ( |50"D with Pi = and s = 0.5, 
1.0, and 2.0. For comparison, we also plot the approximate velocity profile of eq. (|7|). For 
s < ~ 1.5, the flow velocity is larger than the approximate value. This is due to the 1 + fi 2 
dependence of the scattering cross section, which the approximate calculation does not take 
into account. Photons moving at large /i, which carry more momentum along the field line, 
are more likely to scatter, increasing the radiation force and hence the flow velocity. As the 
spectral index increases, there are fewer photons with large energies. Consequently F ra( i is 
smaller and ftp decreases. 

Figure [5| shows the velocity profile for s = 1 and Pi = —1,0, and 1. In perpendicular 
mode, the scattering cross section is independent of /i and the velocity profile is the same as 
the approximate one. As the fraction of parallel mode photons increases, photons at small 
angles to the field are more likely to scatter and the flow velocity increases. 



4.2. Self-Consistent Monte Carlo Calculation 

We calculate self-consistent velocity profiles by iterating calculations of the radiation 
transfer and the hydrodynamics. The procedure is conceptually similar to the one Wang 
and Frank (1981) use to determine the velocity profile and radiation energy density in 
the accretion column of an accretion-powered pulsar. For the radiation transfer, we have 
adapted a Monte Carlo code written for static, slab line-forming regions ( [Wang, Wasserman. 



|fc Salpeter 1988| ; |Wang et al. 1989] ; [Isenbcrg, Lamb, fc Wang 1997] ) to the geometry 
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and physical conditions in an outflow. The code calculates the radiation force, which we 
integrate to determine the velocity profile for the next iteration. We repeat the process 
until the radiation force is zero everywhere, within the statistical uncertainty of the Monte 
Carlo code. 

In addition to the physical processes included in the analytic model of MT, the Monte 
Carlo calculation permits multiple scattering, higher harmonics, and a finite natural line 
width. The analytic model assumes a disk photon source with electrons on the cylinder 
axis, or equivalently, a point photon source with electrons anywhere in the outflow. The 
Monte Carlo code permits electrons anywhere in the outflow and either a point or a disk 
photon source. 

Figure |6| illustrates the effect of these enhancements on the velocity profile. The figure 
shows the velocity profile for a 1/E photon spectrum injected into a line-forming region 
with Th = 0.1 i?*, B Q) i2 = 1.7, and n eo = 10 17 cm -3 . We use polarization averaged cross 
sections in all simulations, except where otherwise indicated. As the figure shows, for a 
Monte Carlo calculation with first harmonic scattering only, a point photon source, and 
T r = 0, the self-consistent velocity profile is very similar to the profile calculated from MT's 
analytic model. The two calculations agree within about 10% at the stellar surface. To 
understand the difference, consider the optical depth of a region which is small enough that 
the magnetic field and flow velocity do not change significantly. The fraction of photons 
with energy E and orientation /i which scatter in this region is f s ~ 1 — exp[— r(E, (/,)]. 
Assuming that first harmonic scattering dominates the radiation transfer, expanding in r, 
and averaging over energy yields, in the frame of reference moving with the flow, 

ff ~ 1 _ p-3/4 n(l+A»°) 
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(57) 



where n is the energy and angle averaged optical depth. When the optical depth is very 
small (n <C 1), is proportional to the cross section, consistent with the single scattering 
approximation used by MT. Since F rac i oc r±, and F ra( i = at the flow velocity, the velocity 
profile is not sensitive to T\ at very small optical depths. However, as the optical depth 



traveling at large angles to the field are more likely to scatter, reducing (3p. 

The finite natural line width reduces the flow velocity even further. We attribute this 
to scattering of photons in the line wings. Wing scattering has the greatest impact at /x = 
where there is no thermal broadening of the cyclotron line. Consequently, the finite natural 
line width increases the number of scatters of photons at large angles to the magnetic field, 
decreasing (3f- 

When photons are injected semi-isotropically at the origin of the hot spot, the mean 
distance from the injection point to the walls of the cylinder is 7r/2 « 1.57 r^. When they 
are injected uniformly across the hot spot, this distance falls slightly, to 4/3 r h rs 1.33 r h , 
resulting in a proportional drop in mean optical depth along the line of sight. Since the 
dependence of the velocity profile on optical depth is second order and the change in optical 
depth is small, the method of injection does not have a significant impact on the profile, as 
Figure |6] shows. 

When higher harmonic scattering is added to the calculation, the flow velocity drops 
further. This is clear from the (1 — /i 2 )™" 1 dependence of the scattering cross sections. 
Photons with small /x are more likely to scatter at higher harmonics, lowering the radiation 
force. 




become significant. Consequently, photons 
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We illustrate the effect of the optical depth on the velocity profile in greater detail in 
Figure |7|. The figure shows velocity profiles for a line-forming region with r% = 0.05 R* and 
-B ,i2 = 1.7. Continuum and first three harmonic scattering with finite natural line width 
are included in the calculation. The injected photon spectrum is the best fit of a double 
power law to the continuum spectrum of GB880205, n x (E) oc E~ s , where s — a± — 0.8479 
for E < E break = 107.8 keV and s = a 2 = 1.199 for E > E break flFreeman 1993| ; see also 
Fenimore et al. 1988] ; |Wang et al. 1989] ). As n eo increases from 10 16 cm~ 3 to 10 17 cm -3 , the 
flow velocity decreases, as expected from eq. (|5"?D . However, as the plasma density increases 
further, to 10 18 cm -3 , the trend reverses. This is because at the larger plasma density there 
is significant angular redistribution of the photons; eq. fl5"7|) does not take this redistribution 
into account. Because of the motion of the plasma, many photons are scattered to large 
values of /i (see Figures || and [15], below). Further scattering of these photons increases F ra d 
and Pf- 



5. Monte Carlo Spectra 

Figure [8] shows the emerging photon flux N(E,fi), in arbitrary units, from a Monte 
Carlo calculation with first harmonic scattering only and zero natural line width. The 
injected photon flux Ni(E,fi) is a power law with s = 1. Injection is at the origin of a hot 
spot with radius = 0.1 i?*. At the stellar surface the magnetic field strength is -B ,i2 — 1-7 
and the electron density is n eo = 10 17 cm -3 . The calculation uses the approximate velocity 
profile of eq. (0). For comparison, we also plot a pure absorption spectrum, 



N abs (E,n) = Ni(E,fj) exp 
where in this case the sum stops at n — 1. 



(58) 
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The motion of the plasma has a significant impact on the spectrum. The line centers 
are shifted from their rest frame values E pa uEb, according to eq. fl3"0|). In addition, the 
high plasma velocity beams scattered photons along the magnetic field. Consequently, when 
viewed perpendicular to the field [panel (a)], the scattered spectrum is nearly identical to 
the absorption spectrum, since photons absorbed at large angles to the field are re-emitted 
at smaller ones. As \x approaches (3p = 1/2, the cutoff energy E c is near the line center 
and plays an important role in determining the properties of the line, as we note in Section 
|3.4j . Since E c falls with altitude, photons with E < E c can escape the line core by diffusing 
upwards until E > E c . The photons then escape the line-forming region without further 
scattering. The escape of photons in the blue wing is responsible for the spike in panel (b). 
This feature becomes broader with increasing /z, due to the enhanced Doppler broadening 
of the line when viewed along the field, and to the large number of scattered photons that 
emerge at high //. As /i approaches unity, the line is almost entirely filled in by scattered 
photons, forming a broad emission-like feature. 

The features in Figure |8| are narrow. The equivalent width of the line in panel (a) 
is Wei = 5.3 keV, which is comparable to the first harmonic line width in the Ginga 
observations of GB880205, Wei = 3.7 keV. The reason the line is so narrow, in spite of the 
variation of the magnetic field with altitude, is that the vast majority of scatters occurs in 
a region close to the surface where the variation of the field is small. We illustrate this in 
Figure |9], which shows the number of scatters as a function of altitude dN s /dz, in arbitrary 
units. We plot dN s /dz for the simulation in Figure || as well as for simulations with first 
three harmonic and continuum scattering and both (3p = 0.5 (1 + //„) and a self-consistent 
velocity profile. Adding the higher harmonics increases the number of scatters. In all cases, 
few scatters occur above z — 1. The decrease in the number of scatters with altitude is 
due to the escape of photons through the sides of the cylindrical line-forming region. The 
decrease is analogous to the decrease in the number of of unscattered photons, described by 
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eq. p. 

Figure |1(] shows that even though the natural line width is small, it has a noticeable 
impact on spectra formed in an outflow, and needs to be included in radiation transfer 
calculations. The figure shows the emerging photon flux from simulations with a 1/E 
spectrum injected at the origin of a hot spot with radius = 0.1 i?*. At the surface 
the magnetic field strength is -B 0j i2 = 1.7 and the electron density is n eo = 10 17 cm" 3 . 
The simulations use self-consistent velocity profiles. When the finite natural line width 
is included, the scattering profile is no longer zero above E c . Consequently, photons in 
the blue wing spike at fi « 1/2 do not escape the line- forming region immediately. They 
continue to scatter to larger values of //, depopulating the spike and filling in the lines at 
higher /x to a greater extent than when T r = 0. 

When the photons are injected uniformly across the hot spot instead of at the origin, 
the lines become shallower. This is due to the reduction in the optical depth seen by 
photons that are injected closer to the edges of the hot spot. We illustrate this in Figure 111] 



for a line-forming region with the same parameters as Figure 10 



Adding continuum and second and third harmonic scattering to the simulation of Figure 



1 1 yields the spectra shown in Figure 12. For comparison, we also plot the corresponding 
spectra for first harmonic scattering only and pure absorption. The higher harmonic 
features in the scattering spectrum are almost identical to the absorption features. This 
is because the scattering cross sections favor Raman scattering at the higher harmonics; 
i.e. most photons absorbed at n=2 and 3 are re-emitted as two or three spawned photons 



at lower harmonics QBussard fc Lamb 1982j ). Consequently the scattering profiles suffice 



to explain the properties of the higher harmonic lines. As \x increases the lines become 
wider due to Doppler broadening, but shallower due to the (1 — /i r2 ) n_1 dependence of 
the scattering cross section. The asymmetric shape is characteristic of lines formed in an 
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outflow — as the altitude increases, both the cyclotron energy and the plasma density 
decrease, resulting in an extended tail in the red wing. This is apparent, for example, in 
the second harmonic line in panel (b). At the first harmonic, the absorption-like features at 
small fj, are shallower and the emission-like features at large \x are more sharply peaked than 
in the simulations with first harmonic scattering only. This is because additional photons 
are present at the first harmonic due to spawning. These photons fill in the absorption-like 
features and build up the emission-like ones. 

In Figure [H| we compare the spectrum of Figure [l2|, calculated from a self-consistent 
velocity profile, to a spectrum based on the approximate velocity profile of eq. (]7|). The 
flow velocity at the surface is (3p = 0.26 in the self-consistent calculation, compared with 
0.50 in the approximation. Because of the smaller plasma velocity near the surface, the line 
centers are shifted closer to their rest frame values. In addition, the lines are considerably 
narrower in the self-consistent calculation. For < fi < 1/8 [panel (a)], Wei = 3.7 keV 
in the self-consistent calculation, compared to 5.1 keV in the approximation. The reason 
for the narrower line is clear from eq. ([27]), which shows that the plasma density n e (z) 
is proportional to the ratio Pfo/ Pf{z)- Figure ^| shows that this ratio is much smaller in 
self-consistent calculations. Consequently, for a given plasma density at z = 0, the plasma 
density at z > is smaller in the self-consistent calculation, as is the optical depth. 



Figure 14 compares the self-consistent calculation of Figure 13, which has temperature 



k B T e = 0.25 E B , to a calculation with k B T e = 0.05 E B . As we discuss in Sections 1372] and 



3.5, the effect of the temperature on the spectra is negligible for the parameters used here. 



Figure [15] illustrates the effect of the surface electron density on the spectra. The 
figure shows the photon flux emerging from a line-forming region with = 0.05 R* and 
S o l2 =1.7. The continuum spectrum of GB880205 is injected uniformly across the hot spot. 
The calculations include scattering in the continuum and the first three harmonics, finite 
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natural line width, and a self-consistent velocity profile. The surface electron densities are 
n eo = 10 16 , 10 17 , and 10 18 cm" 3 . As expected, both the full-width half-maxima and \We\ 
increase with the plasma density. However, even at n eo = 10 18 cm -3 , the equivalent widths 
are modest, with Wei = 5.2 and W E 2 = 7.9 keV for < /j, < 1/8. 

We plot spectra from self-consistent calculations with polarized cross sections in Figure 
16 . The line-forming region has = 0.1 R#, -B D .i2 = 1.7, and n eo = 10 17 cm -3 . The 



calculations include continuum and first three harmonic scattering with finite natural 
line width. The photons, injected uniformly across the hot spot and with the continuum 
spectrum of GB880205, have initial polarizations Pi = —1, 0, and 1. The photons at 
the first harmonic scatter and switch polarization modes many times before escaping. 
Consequently, the properties of the first harmonic lines are not affected very much by the 
use of polarized cross sections or the choice of Pi. The scattering cross sections [eq. (|4C|)1 
explain the properties of the higher harmonic features. Because of the /i r2 dependence of 
the cross section for parallel mode photons, the second and third harmonic lines are very 
shallow near //" = (i.e. \x = (3p ~ 0.3) when Pi = — 1. The depth of the higher harmonic 
lines increases with the fraction of perpendicular mode photons. 



6. Discussion 

Fenimore et al. (1988) report that the two harmonically spaced lines observed by 
Ginga in GB880205 have equivalent widths Wei = 3.7 keV and W E 2 = 9.1 keV. As we show 
in Section |5], cyclotron scattering in an outflow can easily create lines with comparable 
equivalent widths, due to the escape of photons through the sides of the cylindrical 
line-forming region. For example, the first two harmonic lines in the spectrum in Figure 



Top with n eo = 10 cm and 3/8 < /i < 1/2 have equivalent widths Wei = 2.3 keV and 



Wei = 8.8 keV. Clearly, the interpretation of the observed features as cyclotron lines does 
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not rule out burst sources in the galactic corona, provided the photon source is a small 
fraction of the stellar surface. Thus, the cyclotron line interpretation is consistent with the 
BATSE brightness and sky distributions, which suggest that if bursters are galactic, they 
are in the corona. We are currently fitting the outflow model to the Ginga data in greater 
detail by folding the model spectra through the detector response matrices and calculating 
X 2 (see |Eenimore et al. 1981% |Wang et al. 1989| ; |Briggs 19Wj ). 



In addition to generating cyclotron lines comparable to the absorption-like features 
observed by Ginga, the outflow model predicts the formation of emission-like features. 
Some static models predict double-peaked emission-like features, consisting of line shoulders 
on each side of an absorption-like line at certain viewing angles when the photon source is 
embedded in a plasma which is optically thin in the continuum ( [Freeman et al. 1992fc |Araya 



|fc Harding 1996| ; [Isenberg, Lamb, fc Wang 1997|) . But these are quite distinct from the 
szng/e-peaked features that the outflow model predicts. Mazets et al. (1981, 1982) report 
observations of such features by the Konus instruments. The line candidates observed by 
BATSE ( Briggs et al. 19961) also include emission-like features. These observations should 
be compared to the model as soon as the observed spectra and detector response matrices 
become available. 

All the simulations in the present work use plasma densities that are small 
enough so that the line-forming region is optically thin in the continuum, i.e. 
n e ^ o"^ 1 r^ 1 = 1.5 x 10 19 cm" 3 r^\. But for the physical conditions required for 
galactic corona models, the plasma density could be much larger due to the production 
of electron-positron pairs. The cross section for the two photon pair production process 
77 — > e + e~ is ~ cry, corresponding to a pair production rate 
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^~^ = 440c m -. il2rft -,i( r M_)" 1 , (59) 

where L42 = L/10 42 erg s -1 is the total (x-ray plus gamma-ray) burst luminosity, and 
(E) is the mean photon energy In a strong magnetic field, one photon pair production 
(7 — > e + e~) is also significant. The polarization averaged production rate for this process is 



a 



RyB ~ 0.26 — bsin.6 exp 
Ac 



m P c 



3 E bsinO 
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for Eb sin 8/(m e c 2 ) <C 1, and 
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(61) 



for Eb sin 6/(m e c 2 ) ^> 1, where A^ = h/(m e c) is the Compton wavelength ( [l'sai fc Erber 
1974| ; |Daugherty fc Harding l"983| ; [Burns fc Harding 1984fc |Meszaros 199p . Clearly the 



optical depths for both the one and two photon processes are much larger than unity for the 
magnetic fields and photons densities of interest in galactic corona gamma-ray burst models 
( [Schmidt 1978] ; |Brainerd fc Lamb 1987|) . In addition to providing large plasma densities, 
pair production could truncate the spectrum at the pair production threshold, 1 MeV. 
Truncation is not consistent with observations by BATSE ( Band et al. 1993| ), COMPTEL 
( |W inkier et al. 1993| ), and EGRET flSchneid et al. 1992+ |Kwok et al. P^ ; |Sommer et aT 



|1994 ) of photon energies up to 1 GeV with no evidence of a spectral cutoff or rollover at 
1 MeV. 
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The properties of lines formed in an outflow that is optically thick in the continuum 
are not very well known. The full-width half-maxima of the lines increase with optical 
depth like (lnr) 1 / 2 . This line broadening, though logarithmic, is significant if the plasma 
density increases by a factor ^ 100; it might be sufficient to prevent the formation of narrow 
lines. However, parallel mode photons with jj, sa (3p escape the line-forming region easily 
since their resonant scattering cross section is oc /i r2 . Therefore, the spectrum may possess 
narrow cyclotron features when viewed along this line of sight. The photon flux along this 
line of sight is ~ one-tenth the flux at \i — 1 because nearly all photons scatter at these 
optical depths and are beamed along the field. Consequently, if the line of sight is at a large 
angle to the field, the observed flux is reduced and narrow lines are present. If the line 
of sight is at a small angle to the field, the observed flux is large. The absence of narrow 
features in gamma-ray burst spectra calculated from the "CUSP" model, in which the 
spectrum is formed by the Cyclotron Up Scattering Process flVitello fc Dermer 199 1| ; |Coppi 



fc Lamb 1992|) suggests that no narrow cyclotron feature is visible when the line of sight is 
at a small angle to the field. In short, if the cyclotron lines are formed in an outflow that is 
optically thick in the continuum, the presence of narrow lines is inversely correlated with 
the brightness of the burst. This correlation may be masked by selection effects: lines are 
easier to detect in brighter bursts. Whether spectra formed under these circumstances are 
consistent with observations is an open question. Additional simulations are required, using 
cross sections that are appropriate for plasmas that are optically thick in the line wings. 

The plasma densities in corona models may be smaller than eqs. (|59"|)-(|6"1|) suggest (for 
a review see [Higdon fc Lingenfelter 1990|) . One often cited possibility is that the entire 



photon spectrum is produced at the magnetic polar cap, but the photons are beamed into 
a small angle. Another is that the x-rays and gamma-rays are produced in two separate 
physical regions, by two different mechanisms. 
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If a plasma with bulk motion corresponding to 7^ produces the spectrum, then the 
photons are beamed into an angle ~ 7^ . The kinematics of two photon pair production 
requires E\E2{1 — cos 6u) > 2mlc 4 , where E\ and Ei are the energies of the two photons 
and Qyi is the angle between them. A beaming angle ~ 0.1 (i.e. IV ~ 10) raises the 
pair production threshold enough to be consistent with the observed spectra above 1 
MeV ( Schmidt 197$ [Baring h Harding 199$ |fiarding h Baring 1994| |Harding 1994| ). 
Unfortunately, it is unlikely that the observed cyclotron features are produced in a plasma 
with such a large Lorentz factor. According to eq. (|30|), the creation of a cyclotron first 
harmonic at 20 keV when 7^ = 10 and 9 = 0.1 requires a magnetic field Byi ~ 0.2, which 
is about an order of magnitude smaller than the fields we use in the present work. Since 
the ratio of the second and first harmonic scattering cross sections <j 2 /&i oc B, the ratio 
We2/Wei is much smaller for B12 ~ 0.2 than in the Ginga observations. To obtain the 
observed spectrum, the lines must be formed in a plasma in a separate location from the 
plasma responsible for beaming the radiation. The plasma in the line-forming region must 
also be moving more slowly than 7^ = 10. 

This leads to the second possibility, that the gamma-ray burst spectrum consists of 
two components produced in two distinct physical regions: a soft (E ^ 1 MeV) component 
produced near the magnetic pole and a hard component (E ^ 1 MeV), probably produced 
at some distance above the surface (|Lamb 1982| ; Katz 1982, 1994). If this is the case, then 



the photons in the line-forming region have energies below the pair production threshold 
and the plasma density can remain small. Since the region where the gamma-rays are 
produced is larger and has a smaller magnetic field than in models where the gamma-rays 
are produced at the polar cap, there is not sufficient pair production to truncate the 
spectrum at 1 MeV. A recent report by Chernenko and Mitrofanov (1995) of evidence for 
two components in the spectrum of GB881024 is intriguing. The soft component dominates 
the spectrum for E < 250 keV; the hard component is dominant at higher energies. 
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However, the question of exactly where and how the high energy component is produced 
remains open. 

The cyclotron line interpretation of the observed spectral features strongly supports 
a galactic origin for at least some gamma-ray bursts. Recent reports of fading x-ray and 
optical counterparts apparently associated with GB970228 raise hopes that we will soon 
know the distance scale to the burst sources. If it turns out that some bursts are galactic, 
then our understanding of how cyclotron lines are formed at low energies will play an 
important role in constraining models of how the spectrum is formed at all energies. 

We wish to thank Jeff Benensohn, Tomek Bulik, Peter Freeman, Carlo Graziani, 
Cole Miller, Guy Miller, Lucia Muhoz-Franco, Jean Quashnock, Paul Ricker, Jane Wang, 
Rob Wickham, Tom Witten, and Shanqun Zhan for invaluable discussions. We are 
grateful to Tsvi Piran and Stan Woosley for pointing out the necessity of considering the 
impact of pair production on our model. We also wish to thank Corbin Covault, Mark 
Oreglia, and Jonathan Rosner for their careful reading of the manuscript and their helpful 
suggestions. This work was supported in part by NASA Grants NAGW-830, NAGW-1284, 
and NAG5-2868. MI is happy to acknowledge the support of a NASA Graduate Student 
Researchers Fellowship under NASA Grant NGT-5-18. This work partially satisfies the 
requirements for the degree of Ph.D. from the University of Chicago for MI. 



-43- 

REFERENCES 

Adler, S. L. 1971, Ann. Phys., 67, 599 
Alexander, S. G., & Meszaros, P. 1989, ApJ, 344, LI 
Araya, R. A., & Harding, A. K. 1996, ApJ, 463, 33 
Band, D., et al. 1993, ApJ, 413, 281 
Barat, C. 1993, A&AS, 97, 43 

Baring, M. G., & Harding, A. K. 1993, in Proc. 23rd ICRC, Vol. 1, 53 

Bonnell, J. T., & Klebesadel, R. W. 1996, in AIP Conf. Proc. No. 384, Gamma- Ray Bursts, 
3rd Huntsville Symposium, ed. C. Kouveliotou, M. S. Briggs, & G. J. Fishman 
(Woodbury, NY: AIP Press), 977 

Brainerd, J. J., k Lamb, D. Q. 1987, ApJ, 313, 231 

Briggs, M. S. 1996, in AIP Conf. Proc. No. 384, Gamma-Ray Bursts, 3rd Huntsville 

Symposium, ed. C. Kouveliotou, M. S. Briggs, & G. J. Fishman (Woodbury, NY: 
AIP Press), 133 

Briggs, M. S., et al. 1996, in AIP Conf. Proc. No. 384, Gamma- Ray Bursts, 3rd Huntsville 
Symposium, ed. C. Kouveliotou, M. S. Briggs, & G. J. Fishman (Woodbury, NY: 
AIP Press), 153 

Burns, M. L., & Harding, A. K. 1984, ApJ, 285, 747 

Bussard, R. W., & Lamb, F. K. 1982, in AIP Conf. Proc. No. 77, Gamma-Ray Transients 
and Related Astrophysical Phenomena, ed. R. E. Lingenfelter, H. S. Hudson, & 
Diana M. Worrall (New York: American Institute of Physics), 189 



-44- 



Canuto, V., Lodenquai, J., & Ruderman, M. A. 1971, Phys. Rev. D., 3, 2303 
Chanmugam, G. 1992, ARA&A, 30, 143 

Chernenko, A., & Mitrofanov, I. 1993, in Isolated Pulsars, ed. K. A. Van Riper, R. Epstein, 
& C. Ho (Cambridge: Cambridge University Press), 215 

Chernenko, A., & Mitrofanov, I. 1995, MNRAS, 274, 361 

Clark, G. W., et al. 1990, ApJ, 353, 274 

Coppi, P. S., & Lamb, D. Q. 1992, in AIP Conf. Proc. No. 265, Gamma-Ray Bursts, 
Huntsville AL 1991, ed. W. S. Paciesas & G. J. Fishman (New York: American 
Institute of Physics), 257 

Costa, E., et al. 1997, IAU Circ, 6576 

Daugherty, J. K., & Harding, A. K. 1983, ApJ, 273, 761 

Daugherty, J. K., & Ventura, J. 1978, Phys. Rev. D, 18, 1053 

de Haas, W. J., & van Alphen, P. M. 1930, Leiden Comm., 208d 

Fenimore, E. E., et al. 1988, ApJ, 335, L71 

Freeman, P. E. 1993, Private communication 

Freeman, P. E., et al. 1992, in AIP Conf. Proc. No. 265, Gamma-Ray Bursts, Huntsville 
AL 1991, ed. W. S. Paciesas & G. J. Fishman (New York: American Institute of 
Physics), 216 

Freeman, P. E., et al. 1996, in AIP Conf. Proc. No. 366, High Velocity Neutron Stars and 
Gamma-Ray Bursts, ed. R. E. Rothschild & R. E. Lingenfelter (New York: AIP 
Press), 211 



-45 - 



Graziani, C. 1993, ApJ, 412, 351 

Graziani, C, et al. 1992, in Gamma-ray Bursts: Observations, Analyses, and Theories, ed. 
C. Ho, R. I. Epstein, & E. E. Fenimore (Cambridge: Cambridge University Press), 
407 

Groot, P. J., et al. 1997a, IAU Circ., 6584 
Groot, P. J., et al. 1997b, IAU Circ., 6588 
Grove, J. E., et al. 1995, ApJ, 438, L25 

Harding, A. K. 1994, in AIP Conf. Proc. No. 304, The Second Compton Symposium, ed. C. 
Fichtel, N. Gehrels, & J. P. Norris (New York: AIP), 30 

Harding, A. K., & Baring, M. G. 1994, in AIP Conf. Proc. No. 307, Gamma-Ray Bursts, 
Second Workshop, ed. G. J. Fishman, J. J. Brainerd, & K. Hurley (New York: AIP 
Press), 520 

Harding, A. K., k Daugherty, J. K. 1991, ApJ, 374, 687 
Harding, A. K., & Preece, R. D. 1989, ApJ, 338, L21 
Herold, H. 1979, Phys. Rev. D, 19, 2868 
Higdon, J. C, & Lingenfelter, R. E. 1990, ARA&A, 28, 401 

Hueter, G. J. 1984, in AIP Conf. Proc. No. 115, High Energy Transients in Astrophysics, 
ed. S. E. Woosley (New York: American Institute of Physics), 373 

Hueter, G. J. 1987, Ph.D. thesis, University of California, San Diego 

Isenberg, M., Lamb, D. Q., & Wang, J. C. L. 1996, in AIP Conf. Proc. No. 384, 

Gamma-Ray Bursts, 3rd Huntsville Symposium, ed. C. Kouveliotou, M. S. Briggs, 
& G. J. Fishman (Woodbury, NY: AIP Press), 187 



Isenberg, M., Lamb, D. Q., & Wang, J. C. L. 1997, ApJ, submitted 
Katz, J. I. 1982, ApJ, 260, 371 
Katz, J. I. 1994, ApJ, 422, 248 

Klebesadel, R. W., Strong, I. B., & Olson, R. A. 1973, ApJ, 182, L85 

Kwok, P. W., et al. 1993, in AIP Conf. Proc. No. 280, Compton Gamma-ray Observatory, 
ed. M. Friedlander, N. Gehrels, & D. J. Macomb (New York: AIP Press), 855 

Lamb, D. Q. 1982, Ann. N.Y. Acad. Sci., 422, 237 

Lamb, D. Q. 1995, PASP, 107, 1152 

Lamb, D. Q., et al. 1989, in Proc. 14th Texas Symp. on Relativistic Astrophysics, ed. 
E. Fenyves, Vol. 571 (New York: Ann. NY Acad Sci.), 460 

Lamb, D. Q., Wang, J. C. L., k Wasserman, I. 1990, ApJ, 363, 670 

Landau, L. D. 1930, ZPhys, 64, 629 

Landau, L. D., & Lifshitz, E. M. 1975, The Classical Theory of Fields (Oxford: Pergamon 
Press) 

Landau, L. D., & Lifshitz, E. M. 1980, Statistical Physics, Part I (Oxford: Pergamon Press) 
Langer, S. 1981, Phys. Rev. D., 23, 328 

Makishima, K., & Mihara, T. 1992, in Frontiers of X-Ray Astronomy, ed. Y. Tanaka & K. 
Koyama (Tokyo: Universal Academy Press), 23 

Mazets, E. P., et al. 1981, Nature, 290, 378 

Mazets, E. P., et al. 1982, Ap&SS, 82, 261 



-47- 



Mazets, E. P., et al. 1983, in AIP Conf. Proc. No. 101, Positron-electron pairs in 

astrophysics, ed. M. L. Burns, A. K. Harding, & R. Ramaty (New York: AIP Press), 
36 

Mazets, E. P., et al. 1996, in AIP Conf. Proc. No. 384, Gamma- Ray Bursts, 3rd Huntsville 
Symposium, ed. C. Kouveliotou, M. S. Briggs, & G. J. Fishman (Woodbury, NY: 
AIP Press), 492 

Meegan, C, et al. 1996, IAU Circ, 6518 

Meszaros, P. 1992, High Energy Radiation from Magnetized Neutron Stars (Chicago: 
University of Chicago Press) 

Metzger, M. R., et al. 1997a, IAU Circ, 6588 

Metzger, M. R., et al. 1997b, IAU Circ, 6631 

Mihara, T. 1995, Ph.D. thesis, University of Tokyo 

Miller, G. S., et al. 1991, Phys. Rev. Lett., 66, 1395 

Miller, G. S., et al. 1992, in Gamma-ray Bursts: Observations, Analyses, and Theories, ed. 
C. Ho, R. I. Epstein, & E. E. Fenimore (Cambridge: Cambridge University Press), 
215 

Mitrofanov, I. G., & Tsygan, A. I. 1982, Ap&SS, 84, 35 
Murakami, T., et al. 1988, Nature, 335, 234 

Nelson, R. W., Salpeter, E. E., & Wasserman, I. 1993, ApJ, 418, 874 

Nishimura, O. 1994, PASJ, 46, 45 

Nishimura, O., & Ebisuzaki, T. 1992, PASJ, 44, 109 



-48- 

Paczynski, B. 1995, PASP, 107, 1167 

Rybicki, G. B., & Lightman, A. P. 1979, Radiative Processes in Astrophysics (New York: 
John Wiley and Sons) 

Sahu, K., et al. 1997, IAU Circ, 6606 

Schmidt, W. K. H. 1978, Nature, 271, 525 

Schneid, E., et al. 1992, A&A, 255, L13 

Sommer, M., et al. 1994, ApJ, 422, L63 

Spitzer, L. 1962, Physics of Fully Ionized Gases (New York: Interscience Publishers) 

Ternov, V. G., Begrov, V. G., & Zhukovskii, V. Ch. 1966, Moscow Univ. Bull., 21, 21 

Triimper, J., et al. 1978, ApJ, 219, L105 

Tsai, W., & Erber, T. 1974, Phys. Rev. D, 10, 492 

Tsygan, A. I. 1981, Ap&SS, 77, 187 

Ventura, J. 1979, Phys. Rev. D, 19, 1684 

Visvanathan, N., & Wickramasinghe, D. T. 1979, Nature, 281, 47 

Vitello, P., & Dermer, C. D. 1991, ApJ, 374, 668 

Wang, J. C. L., et al. 1989, Phys. Rev. Lett., 63, 1550 

Wang, J. C. L., Wasserman, I., & Lamb, D. Q. 1993, ApJ, 414, 815 

Wang, J. C. L., Wasserman, I. M., & Salpeter, E. E. 1988, ApJS, 68, 735 

Wang, Y.-M., & Frank, J. 1981, A&A, 92, 255 



-49- 



Wasserman, L, k Salpeter, E. 1980, ApJ, 241, 1107 
Wheaton, W. A., et al. 1979, Nature, 282, 240 

Winkler, C, et al. 1993, in AIP Conf. Proc. No. 280, Compton Gamma-ray Observatory, 
ed. M. Friedlander, N. Gehrels, & D.J. Macomb (New York: AIP Press), 845 

Yoshida, A., et al. 1992, in Gamma-ray Bursts: Observations, Analyses, and Theories, ed. 
C. Ho, R. I. Epstein, & E. E. Fenimore (Cambridge: Cambridge University Press), 
399 

Zheleznyakov, V. V., & Serber, A. V. 1994, Sp. Sci. Rev., 68, 275 
Zheleznyakov, V. V., & Serber, A. V. 1995, Adv. Sp. Res., 16, 77 



This manuscript was prepared with the AAS IATjtN. macros v4.0. 

Fig. 1. — Geometry of the Line-Forming Region. Photons are injected at a hot spot with 
radius r h located at the magnetic pole of a neutron star with dipole field B. They scatter 
with the electrons in the line-forming region, accelerating them along the field lines. (5 is the 
velocity of an individual electron and (3f{z) is the bulk velocity of the flow, where z is the 
altitude above the stellar surface. 9 denotes the angle between the photon direction and the 
magnetic field direction. The angular size of the hot spot is o (z) 
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Fig. 2. — First harmonic resonant photon energy E\ as a function of electron momentum 
p for B12 = 1.7 and several photon orientations /i. For a given [i and E\ the resonant 
condition has two roots p^ (solid lines) and p + (dotted lines). Note the cutoff energy E c 
where p_ = p + and (3 = fx (i.e. /i r = 0). Physical roots do not exist for E > E c . The values 
of p corresponding to (5 — and j3 = 0.5 are indicated by dashed lines. 

Fig. 3. — Radiation-driven acceleration of an electron-proton pair in a strong magnetic field, 
according to the model of Mitrofanov and Tsygan (1982). The particles are injected at the 
origin of a hot spot with radius r^. Radiation with a 1/E spectrum and X-ray luminosity 
L x = 10 40 erg s _1 is injected semi-isotropically and uniformly across the hot spot. Left: 
velocity (3 as a function of altitude z. Right: Lorentz factor 7. The solid lines are the (5 and 
7 of a pair with zero initial velocity. The dotted lines indicate the values of f3 and 7 where 
the radiation force F rad = 0. 

Fig. 4. — Effect of the spectral index s on the velocity profile. Flow velocity as a function of 
altitude calculated from the model of Mitrofanov and Tsygan (1982) for s=0.5 (dotted line), 
1.0 (dashed line), and 2.0 (dot-dash) and unpolarized scattering. The approximate profile 
Pf = 0.5 (1 + Ho) is shown for comparison (solid line). 

Fig. 5. — Effect of the photon polarization on the velocity profile. Flow velocity as a function 
of altitude calculated from the model of Mitrofanov and Tsygan (1982) for spectral index 
s=l and polarized scattering. The fraction of injected photons in parallel mode is Pi = —1.0 
(dashed line), 0.0 (dotted line), and 1.0 (solid line). 
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Fig. 6. — Self-consistent Monte Carlo calculations of velocity profiles for a 1/E spectrum 
injected into a line- forming region with hot spot radius = 0.1 dipole field B 0ti2 = 1.7, 
and surface electron density n efi = 10 17 cm -3 . The profile is shown for first harmonic 
scattering, point photon source, and natural line width r r =0 (stars); first harmonic 
scattering, point source, and T r =finite (squares); first harmonic scattering, disk source, 
and r r =finite (dots); and continuum and first three harmonic scattering, disk source, and 
r r =finite (triangles). The profile calculated from the model of Mitrofanov and Tsygan (1982) 
is shown for comparison (solid line). 

Fig. 7. — Effect of electron density on the velocity. Velocity profiles for a GB880205 
continuum spectrum injected uniformly across a hot spot with radius r h = 0.05 The 
line-forming region has dipole field B yi = 1-7, electron temperature ksTl = 0.25 Eb, a 
self-consistent velocity profile, unpolarized continuum and first three harmonic scattering 
with finite natural line width, and surface electron densities n e>0 = 10 16 cm -3 (squares), 
n et0 = 10 17 cm -3 (dots), and n £t0 = 10 18 cm" 3 (triangles). The profile calculated from the 
model of Mitrofanov and Tsygan (1982) is shown for comparison (solid line). 

Fig. 8. — Monte Carlo emerging photon number spectra for several viewing angles for a 
1/E continuum spectrum injected at the origin of a hot spot with radius r h = 0.1 R*. The 
line-forming region has a dipole field B Q:12 = 1.7, electron temperature /ceT e f = 0.25 E B , 
surface electron density n eiG = 10 17 cm -3 , velocity profile f3p = 0.5 (1 + ji ), and unpolarized 
first harmonic scattering with zero natural line width (solid lines). Pure absorption spectra 
(dashed lines) are shown for comparison. 
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Fig. 9. — Number of scatters per photon per unit distance dN s /dz, in arbitrary units, as a 
function of altitude above the stellar surface z/rh for a 1/E continuum spectrum injected into 
a line-forming region with dipole field -B ,i2 = 1-7, electron temperature ksTl = 0.25 E B , 
surface electron density n eiG = 10 17 cm -3 , and unpolarized scattering. Calculations are 
shown for continuum and first three harmonic scattering with self-consistent (solid line) 
and /3p = 0.5 (1 + /i ) (dotted line) velocity profiles, and first harmonic scattering with a 
Pf = 0.5 (1 + [/, ) (dashed line) velocity profile. 



Fig. 10. — Effect of finite natural line width. Monte Carlo emerging photon number spectra 
for several viewing angles for a 1/E continuum spectrum injected at the origin of a hot 
spot with radius r\ = 0.1 R*. The line-forming region has dipole field B ^ 2 = 1-7, electron 
temperature h B T* = 0.25 E B , surface electron density n ej0 = 10 17 cm -3 , a self-consistent 
velocity profile, and unpolarized first harmonic scattering with finite (solid lines) and zero 
(dotted lines) natural line width. 



Fig. 11. — Effect of disk source. Monte Carlo emerging photon number spectra for several 
viewing angles for a 1/E continuum spectrum injected uniformly across a hot spot with 
radius r\ = 0.1 i?* (solid lines) and at the origin (dotted lines). The line-forming region 
has dipole field B Q:12 = 1.7, electron temperature k B Tl = 0.25 E B , surface electron density 
n 6i0 = 10 17 cm~ 3 , a self-consistent velocity profile, and unpolarized first harmonic scattering 
with finite natural line width. 
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Fig. 12. — Effect of higher harmonic scattering. Monte Carlo emerging photon number 
spectra for several viewing angles for a 1/E continuum spectrum injected uniformly across 
a hot spot with radius r\ = 0.1 i?*. The line-forming region has dipole field -B ,i2 = 1-7, 
electron temperature fceT e f = 0.25 E B , surface electron density n £j0 = 10 17 cm -3 , a self- 
consistent velocity profile, and unpolarized continuum and first three harmonic scattering 
with finite natural line width (solid lines). Spectra for first harmonic scattering (dotted lines) 
and pure absorption (dashed lines) are shown for comparison. 



Fig. 13. — Monte Carlo emerging photon number spectra for several viewing angles for a 1/E 
continuum spectrum injected uniformly across a hot spot with r h = 0.1 The line-forming 
region has dipole field B ^ 2 = 1-7, electron temperature ksT^ = 0.25 Eb, surface electron 
density n efi = 10 17 cm -3 , unpolarized continuum and first three harmonic scattering with 
finite natural line width, and a self-consistent (solid lines) and (3p = 0.5 (1 + fi Q ) (dotted 
lines) velocity profile. 



Fig. 14. — Effect of electron temperature. Monte Carlo emerging photon number spectra 
for several viewing angles for a 1/E continuum spectrum injected uniformly across a hot 
spot with radius r h = 0.1 i?*. The line-forming region has dipole field B ol2 = 1.7, surface 
electron density n ej0 = 10 17 cm~ 3 , a self-consistent velocity profile, unpolarized continuum 
and first three harmonic scattering with finite natural line width, and electron temperatures 
k B Tl = 0.05 E B (solid lines) and 0.25 E B (dotted lines). 
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Fig. 15. — Effect of electron density. Monte Carlo emerging photon number spectra for 
several viewing angles for a GB880205 continuum spectrum injected uniformly across a hot 
spot with radius r h = 0.05 R*. The line-forming region has dipole field B o l2 = 1.7, electron 
temperature /csT e f = 0.25 Eb, a self-consistent velocity profile, unpolarized continuum and 
first three harmonic scattering with finite natural line width, and surface electron densities 
n et0 = 10 16 cm" 3 (solid lines), n efi = 10 17 cm -3 (dotted lines), and n ej0 = 10 18 cm -3 (dashed 
lines). 



Fig. 16. — Effect of photon polarization. Monte Carlo emerging photon number spectra 
for several viewing angles for a GB880205 continuum spectrum injected uniformly across 
a hot spot with radius r h = 0.1 R*. The line-forming region has dipole field B 0:12 = 1.7, 
electron temperature hsT* = 0.25 E B , surface electron density n ej0 = 10 17 cm~ 3 , a self- 
consistent velocity profile, and polarized continuum and first three harmonic scattering with 
finite natural line width. The fraction of injected photons in parallel mode is Pi = —1.0 
(dashed lines), 0.0 (dotted lines), and 1.0 (solid lines). 
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Table 1. Summary of calculations of cyclotron line formation in radiation-driven outflows 

associated with gamma-ray bursts. 





Miller et al. 


Chernenko & 


Present 




1991, 1992 


Mitrofanov 1993 


Work 


Process 


Absorption 


Absorption 


Scattering 


Continuum Scat. 


No 


No 


Yes 


Natural line width 





finite 


finite 


Thermal broadening 


Yes 


No 


Yes 


Polarization 


No 


No 


Yes 


Self-consistency 


No 


No 


Yes 
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